Ago2 and a miRNA reduce Topoisomerase 1 for enhancing DNA cleavage in antibody diversification by activation-induced cytidine deaminase.

Activation-induced cytidine deaminase (AID) is the essential enzyme for imprinting immunological memory through class switch recombination (CSR) and somatic hypermutation (SHM) of the immunoglobulin (Ig) gene. AID-dependent reduction of Topoisomerase 1 (Top1) promotes DNA cleavage that occurs upon Ig gene diversification, whereas the mechanism behind AID-induced Top1 reduction remains unclear. Here, we clarified the contribution of the microRNA-Ago2 complex in AID-dependent Top1 decrease. Ago2 binds to Top1 3'UTR with two regions of AID-dependent Ago2-binding sites (5'- and 3'dABs). Top1 3'UTR knockout (3'UTRKO) in B lymphoma cells leads to decreases in DNA break efficiency in the IgH gene accompanied by a reduction in CSR and SHM frequencies. Furthermore, AID-dependent Top1 protein reduction and Ago2-binding to Top1 mRNA are down-regulated in 3'UTRKO cells. Top1 mRNA in the highly translated fractions of the sucrose gradient is decreased in an AID-dependent and Top1 3'UTR-mediated manner, resulting in a decrease in Top1 protein synthesis. Both AID and Ago2 localize in the mRNA-binding protein fractions and they interact with each other. Furthermore, we found some candidate miRNAs which possibly bind to 5'- and 3'dAB in Top1 mRNA. Among them, miR-92a-3p knockdown induces the phenotypes of 3'UTRKO cells to wild-type cells whereas it does not impact on 3'UTRKO cells. Taken together, the Ago2-miR-92a-3p complex will be recruited to Top1 3'UTR in an AID-dependent manner and posttranscriptionally reduces Top1 protein synthesis. These consequences cause the increase in a non-B-DNA structure, enhance DNA cleavage by Top1 in the Ig gene and contribute to immunological memory formation.

De Novo Genome Assembly of Japanese Black Cattle as Model of an Economically Relevant Animal.

A genetic analysis of Japanese Black cattle using short reads and guided by the reference genome from Western breeds would miss the structural variation and/or other unique characteristics of Japanese Black cattle. To overcome this difficulty, a de novo genome assembly independent from the reference genome is required. This chapter describes the technical developments, with respect to both experimental and bioinformatics procedures, including the use of short and long reads, required for de novo genome assembly of Japanese Black cattle.

Comprehensive analysis of 124 transcriptomes from 31 tissues in developing, juvenile, and adult Japanese Black cattle.

Omic analyses of economically important animals, including Japanese Black cattle, are currently underway worldwide. In particular, tissue and developmental stage-specific transcriptome characterization is essential for understanding the molecular mechanisms underlying the phenotypic expression of genetic disorders and economic traits. Here, we conducted a comprehensive analysis of 124 transcriptomes across 31 major tissues from fetuses, juvenile calves, and adult Japanese Black cattle using short-read sequencing. We found that genes exhibiting high tissue-specific expression tended to increase after 60 days from fertilization and significantly reflected tissue-relevant biology. Based on gene expression variation and inflection points during development, we categorized gene expression patterns as stable, increased, decreased, temporary, or complex in each tissue. We also analysed the expression profiles of causative genes (e.g. SLC12A1, ANXA10, and MYH6) for genetic disorders in cattle, revealing disease-relevant expression patterns. In addition, to directly analyse the structure of full-length transcripts without transcript reconstruction, we performed RNA sequencing analysis of 22 tissues using long-read sequencing and identified 232 novel non-RefSeq isoforms. Collectively, our comprehensive transcriptomic analysis can serve as an important resource for the biological and functional interpretation of gene expression and enable the mechanistic interpretation of genetic disorders and economic traits in Japanese Black cattle.

Identification of deleterious recessive haplotypes and candidate deleterious recessive mutations in Japanese Black cattle.

Intensive use of a few elite sires has increased the risk of the manifestation of deleterious recessive traits in cattle. Substantial genotyping data gathered using single-nucleotide polymorphism (SNP) arrays have identified the haplotypes with homozygous deficiency, which may compromise survival. We developed Japanese Black cattle haplotypes (JBHs) using SNP array data (4843 individuals) and identified deleterious recessive haplotypes using exome sequencing of 517 sires. We identified seven JBHs with homozygous deficiency. JBH_10 and JBH_17 were associated with the resuming of estrus after artificial insemination, indicating that these haplotypes carried deleterious mutations affecting embryonic survival. The exome data of 517 Japanese Black sires revealed that AC_000165.1:g.85341291C>G of IARS in JBH_8_2, AC_000174.1:g.74743512G>T of CDC45 in JBH_17, and a copy variation region (CNVR_27) of CLDN16 in JBH_1_1 and JBH_1_2 were the candidate mutations. A novel variant AC_000174.1:g.74743512G>T of CDC45 in JBH_17 was located in a splicing donor site at a distance of 5 bp, affecting pre-mRNA splicing. Mating between heterozygotes of JBH_17 indicated that homozygotes carrying the risk allele died around the blastocyst stage. Analysis of frequency of the CDC45 risk allele revealed that its carriers were widespread throughout the tested Japanese Black cattle population. Our approach can effectively manage the inheritance of recessive risk alleles in a breeding population.

A 44-kb deleted-type copy number variation is associated with decreasing complement component activity and calf mortality in Japanese Black cattle.

BACKGROUND: Calf mortality generally occurs in calves prior to weaning, which is a serious problem in cattle breeding. Several causative variants of monogenic Mendelian disorders in calf mortality have been identified, whereas genetic factors affecting the susceptibility of calves to death are not well known. To identify variants associated with calf mortality in Japanese Black cattle, we evaluated calf mortality as a categorical trait with a threshold model and performed a genome-wide copy number variation (CNV) association study on calf mortality. RESULTS: We identified a 44-kb deleted-type CNV ranging from 103,317,687 to 103,361,802 bp on chromosome 5, which was associated with the mortality of 1-180-day-old calves. The CNV harbored C1RL, a pseudogene, and an IncRNA localized in the C1R and C1S gene cluster, which is a component of the classical complement activation pathway for immune complexes for infectious pathogens. The average complement activity in CNVR_221 homozygotes at postnatal day 7 was significantly lower than that of wild-type animals and heterozygotes. The frequency of the risk allele in dead calves suffering from diarrhea and pneumonia and in healthy cows was 0.35 and 0.28, respectively (odds ratio = 2.2, P = 0.016), suggesting that CNVR_221 was associated with the mortality of Japanese Black calves suffering from an infectious disease. CONCLUSIONS: This study identified a deleted-type CNV associated with the mortality of 1-180-day-old calves. The complement activity in CNVR_221 homozygotes was significantly lower than that in heterozygotes and wild type animals. The frequency of the risk allele was higher in dead calves suffering from an infectious disease than in healthy cows. These results suggest that the existence of CNVR_221 in calves could be attributed to a reduction in complement activity, which in turn leads to susceptibility to infections. Thus, the risk allele could serve as a useful marker to reduce the mortality of infected Japanese Black calves.

New type of encephalomyelitis responsive to trimethoprim/sulfamethoxazole treatment in Japan.

OBJECTIVE: To determine the causative pathogen and investigate the effective treatment of a new type of encephalomyelitis with an unknown pathogen in Japan and report the preliminary ultrastructural and genomic characterization of the causative agent. METHODS: From 2005 to 2012, we treated 4 Japanese patients with geographic clustering and comparable clinical features, serum/CSF cytology, and radiologic findings. Brain biopsy was conducted in all patients to analyze neuropathologic changes by histology, and electron microscopy was applied to reveal the features of the putative pathogen. Genomic DNA was obtained from the affected brain tissues and CSF, and an unbiased high-throughput sequencing approach was used to screen for specific genomic sequences indicative of the pathogen origin. RESULTS: All patients exhibited progressive dementia with involuntary tongue movements. Cytologic examination of CSF revealed elevated mononuclear cells. Abnormal MRI signals were observed in temporal lobes, subcortical white matter, and spinal cord. Biopsied brain tissue exhibited aggregated periodic acid-Schiff-positive macrophages and 2-7 µm diameter round/oval bodies without nuclei or cell walls scattered around the vessels. Unbiased high-throughput sequencing identified more than 100 archaea-specific DNA fragments. All patients were responsive to trimethoprim/sulfamethoxazole (TMP-SMX) plus corticosteroid therapy. CONCLUSIONS: We report 4 cases of encephalomyelitis due to an unknown pathogen. On the basis of ultrastructural and genomic studies, we propose a new disease entity resulting from a causative pathogen having archaeal features. TMP-SMX therapy was effective against this new type of encephalomyelitis.

DB-AT: a 2015 update to the Full-parasites database brings a multitude of new transcriptomic data for apicomplexan parasites.

The previous release of our Full-parasites database (http://fullmal.hgc.jp/) brought enhanced functionality, an expanded full-length cDNA content, and new RNA-Seq datasets from several important apicomplexan parasites. The 2015 update witnesses the major shift in the databases content with focus on diverse transcriptomes of the apicomplexan parasites. The content of the database was substantially enriched with transcriptome information for new apicomplexan parasites. The latest version covers a total of 17 species, with addition of our newly generated RNA-Seq data of a total of 909,150,388 tags. Moreover, we have generated and included two novel and unique datasets, which represent diverse nature of transcriptomes in individual parasites in vivo and in vitro. One is the data collected from 116 Indonesian patients infected with Plasmodium falciparum. The other is a series of transcriptome data collected from a total of 38 single cells of P. falciparum cultured in vitro. We believe that with the recent advances our database becomes an even better resource and a unique platform in the analysis of apicomplexan parasites and their interaction with their hosts. To adequately reflect the recent modifications and the current content we have changed the database name to DB-AT--DataBase of Apicomplexa Transcriptomes.

DBTSS as an integrative platform for transcriptome, epigenome and genome sequence variation data.

DBTSS (http://dbtss.hgc.jp/) was originally constructed as a collection of uniquely determined transcriptional start sites (TSSs) in humans and some other species in 2002. Since then, it has been regularly updated and in recent updates epigenetic information has also been incorporated because such information is useful for characterizing the biological relevance of these TSSs/downstream genes. In the newest release, Release 9, we further integrated public and original single nucleotide variation (SNV) data into our database. For our original data, we generated SNV data from genomic analyses of various cancer types, including 97 lung adenocarcinomas and 57 lung small cell carcinomas from Japanese patients as well as 26 cell lines of lung cancer origin. In addition, we obtained publically available SNV data from other cancer types and germline variations in total of 11,322 individuals. With these updates, users can examine the association between sequence variation pattern in clinical lung cancers with its corresponding TSS-seq, RNA-seq, ChIP-seq and BS-seq data. Consequently, DBTSS is no longer a mere storage site for TSS information but has evolved into an integrative platform of a variety of genome activity data.

Aberrant transcriptional regulations in cancers: genome, transcriptome and epigenome analysis of lung adenocarcinoma cell lines.

Here we conducted an integrative multi-omics analysis to understand how cancers harbor various types of aberrations at the genomic, epigenomic and transcriptional levels. In order to elucidate biological relevance of the aberrations and their mutual relations, we performed whole-genome sequencing, RNA-Seq, bisulfite sequencing and ChIP-Seq of 26 lung adenocarcinoma cell lines. The collected multi-omics data allowed us to associate an average of 536 coding mutations and 13,573 mutations in promoter or enhancer regions with aberrant transcriptional regulations. We detected the 385 splice site mutations and 552 chromosomal rearrangements, representative cases of which were validated to cause aberrant transcripts. Averages of 61, 217, 3687 and 3112 mutations are located in the regulatory regions which showed differential DNA methylation, H3K4me3, H3K4me1 and H3K27ac marks, respectively. We detected distinct patterns of aberrations in transcriptional regulations depending on genes. We found that the irregular histone marks were characteristic to EGFR and CDKN1A, while a large genomic deletion and hyper-DNA methylation were most frequent for CDKN2A. We also used the multi-omics data to classify the cell lines regarding their hallmarks of carcinogenesis. Our datasets should provide a valuable foundation for biological interpretations of interlaced genomic and epigenomic aberrations.

The Babesia bovis gene and promoter model: an update from full-length EST analysis.

BACKGROUND: Babesia bovis is an apicomplexan parasite that causes babesiosis in infected cattle. Genomes of pathogens contain promising information that can facilitate the development of methods for controlling infections. Although the genome of B. bovis is publically available, annotated gene models are not highly reliable prior to experimental validation. Therefore, we validated a preproposed gene model of B. bovis and extended the associated annotations on the basis of experimentally obtained full-length expressed sequence tags (ESTs). RESULTS: From in vitro cultured merozoites, 12,286 clones harboring full-length cDNAs were sequenced from both ends using the Sanger method, and 6,787 full-length cDNAs were assembled. These were then clustered, and a nonredundant referential data set of 2,115 full-length cDNA sequences was constructed. The comparison of the preproposed gene model with our data set identified 310 identical genes, 342 almost identical genes, 1,054 genes with potential structural inconsistencies, and 409 novel genes. The median length of 5' untranslated regions (UTRs) was 152 nt. Subsequently, we identified 4,086 transcription start sites (TSSs) and 2,023 transcriptionally active regions (TARs) by examining 5' ESTs. We identified ATGGGG and CCCCAT sites as consensus motifs in TARs that were distributed around -50 bp from TSSs. In addition, we found ACACA, TGTGT, and TATAT sites, which were distributed periodically around TSSs in cycles of approximately 150 bp. Moreover, related periodical distributions were not observed in mammalian promoter regions. CONCLUSIONS: The observations in this study indicate the utility of integrated bioinformatics and experimental data for improving genome annotations. In particular, full-length cDNAs with one-base resolution for TSSs enabled the identification of consensus motifs in promoter sequences and demonstrated clear distributions of identified motifs. These observations allowed the illustration of a model promoter composition, which supports the differences in transcriptional regulation frameworks between apicomplexan parasites and mammals.

Construction of mate pair full-length cDNAs libraries and characterization of transcriptional start sites and termination sites.

To identify and characterize transcript structures ranging from transcriptional start sites (TSSs) to poly(A)-addition sites (PASs), we constructed and analyzed human TSS/PAS mate pair full-length cDNA libraries from 14 tissue types and four cell lines. The collected information enabled us to define TSS cluster (TSC) and PAS cluster (PAC) relationships for a total of 8530/9400 RefSeq genes, as well as 4251/5618 of their putative alternative promoters/terminators and 4619/4605 intervening transcripts, respectively. Analyses of the putative alternative TSCs and alternative PACs revealed that their selection appeared to be mostly independent, with rare exceptions. In those exceptional cases, pairs of transcript units rarely overlapped one another and were occasionally separated by Rad21/CTCF. We also identified a total of 172 similar cases in which TSCs and PACs spanned adjacent but distinct genes. In these cases, different transcripts may utilize different functional units of a particular gene or of adjacent genes. This approach was also useful for identifying fusion gene transcripts in cancerous cells. Furthermore, we could construct cDNA libraries in which 3'-end mate pairs were distributed randomly over the transcripts. These libraries were useful for assembling the internal structure of previously uncharacterized alternative promoter products, as well as intervening transcripts.

Comparative genome analysis of three eukaryotic parasites with differing abilities to transform leukocytes reveals key mediators of Theileria-induced leukocyte transformation.

We sequenced the genome of Theileria orientalis, a tick-borne apicomplexan protozoan parasite of cattle. The focus of this study was a comparative genome analysis of T. orientalis relative to other highly pathogenic Theileria species, T. parva and T. annulata. T. parva and T. annulata induce transformation of infected cells of lymphocyte or macrophage/monocyte lineages; in contrast, T. orientalis does not induce uncontrolled proliferation of infected leukocytes and multiplies predominantly within infected erythrocytes. While synteny across homologous chromosomes of the three Theileria species was found to be well conserved overall, subtelomeric structures were found to differ substantially, as T. orientalis lacks the large tandemly arrayed subtelomere-encoded variable secreted protein-encoding gene family. Moreover, expansion of particular gene families by gene duplication was found in the genomes of the two transforming Theileria species, most notably, the TashAT/TpHN and Tar/Tpr gene families. Gene families that are present only in T. parva and T. annulata and not in T. orientalis, Babesia bovis, or Plasmodium were also identified. Identification of differences between the genome sequences of Theileria species with different abilities to transform and immortalize bovine leukocytes will provide insight into proteins and mechanisms that have evolved to induce and regulate this process. The T. orientalis genome database is available at http://totdb.czc.hokudai.ac.jp/.

Monitoring endoplasmic reticulum stress responsive mRNAs by RNA sequencing.

Gene expression profile upon endoplasmic reticulum (ER) stress was analyzed by deep shotgun sequencing of mRNAs (DSSR) using RNAs from polysomes or cytoplasm of the HT29 cell. Two time points, 4h after tunicamycin treatment when IRE1α signaling pathway is active and 16h after the treatment when it is inactive, were used. There was a transient decrease in the proportion of shorter mRNA species (<1000bp) in polysome, while it increased transiently in the cytoplasm. Despite such an overall change and decrease in total amount of polysomes, the majority of the 6966 genes analyzed had less than 2 fold change in their expressions. We searched for the genes whose expression was elevated by 2 folds or more in both polysome and cytoplasm and confirmed the results with RT-PCR. There were 7 genes elevated only at 4h (Group I), 20 genes only at 16h (Group II) and 7 genes both at 4 and 16h (Group III). There were 3 genes involved in ribosomal RNA biogenesis in Group I and 2 genes involved mTOR control in Group III. This was consistent with the concept that the ribosome is the essential site for managing ER stress. DSSR is a useful tool for the search of candidates of ER stress responsive genes.

A pilot study on developing mucosal vaccine against alveolar echinococcosis (AE) using recombinant tetraspanin 3: Vaccine efficacy and immunology.

BACKGROUND: We have previously evaluated the vaccine efficacies of seven tetraspanins of Echinococcus multilocularis (Em-TSP1-7) against alveolar echinococcosis (AE) by subcutaneous (s.c.) administration with Freund's adjuvant. Over 85% of liver cyst lesion number reductions (CLNR) were achieved by recombinant Em-TSP1 (rEm-TSP1) and -TSP3 (rEm-TSP3). However, to develop an efficient and safe human vaccine, the efficacy of TSP mucosal vaccines must be thoroughly evaluated. METHODOLOGY/PRINCIPAL FINDINGS: rEm-TSP1 and -TSP3 along with nontoxic CpG ODN (CpG oligodeoxynucleotides) adjuvant were intranasally (i.n.) immunized to BALB/c mice and their vaccine efficacies were evaluated by counting liver CLNR (experiment I). 37.1% (p < 0.05) and 62.1% (p < 0.001) of CLNR were achieved by these two proteins, respectively. To study the protection-associated immune responses induced by rEm-TSP3 via different immunization routes (i.n. administration with CpG or s.c. immunization with Freund's adjuvant), the systemic and mucosal antibody responses were detected by ELISA (experiment II). S.c. and i.n. administration of rEm-TSP3 achieved 81.9% (p < 0.001) and 62.8% (p < 0.01) CLNR in the liver, respectively. Both the immunization routes evoked strong serum IgG, IgG1 and IgG2α responses; i.n. immunization induced significantly higher IgA responses in nasal cavity and intestine compared with s.c. immunization (p < 0.001). Both immunization routes induced extremely strong liver IgA antibody responses (p < 0.001). The Th1 and Th2 cell responses were assessed by examining the IgG1/IgG2α ratio at two and three weeks post-immunization. S.c. immunization resulted in a reduction in the IgG1/IgG2α ratio (Th1 tendency), whereas i.n. immunization caused a shift from Th1 to Th2. Moreover, immunohistochemistry showed that Em-TSP1 and -TSP3 were extensively located on the surface of E. multilocularis cysts, protoscoleces and adult worms with additional expression of Em-TSP3 in the inner part of protoscoleces and oncospheres. CONCLUSIONS: Our study indicated that i.n. administration of rEm-TSP3 with CpG is able to induce both systemic and local immune responses and thus provides significant protection against AE.

Construction and analysis of full-length cDNA library of Cryptosporidium parvum.

A full-length cDNA library was constructed from the sporozoite of Cryptosporidium parvum. Normalized clones were subjected to Solexa shotgun sequencing, and then complete sequences for 1066 clones were reconfigured. Detailed analyses of the sequences revealed that 13.5% of the transcripts were spliced; the average and median 5' UTR lengths were 213.5 and 122 nucleotides, respectively. There were 148 inconsistencies out of 562 examined genes between the experimentally described cDNA sequence and the predicted sequence from its genome. In addition, we identified 118 sequences that had little homology against annotated genes of C. parvum as prospective candidates for addable genes. These observations should improve the reliability of C. parvum transcriptome and provide a versatile resource for further studies.

Full-parasites: database of full-length cDNAs of apicomplexa parasites, 2010 update.

Full-Parasites (http://fullmal.hgc.jp/) is a transcriptome database of apicomplexa parasites, which include Plasmodium and Toxoplasma species. The latest version of Full-Parasites contains a total of 105,786 EST sequences from 12 parasites, of which 5925 full-length cDNAs have been completely sequenced. Full-Parasites also contain more than 30 million transcription start sites (TSS) for Plasmodium falciparum (Pf) and Toxoplasma gondii (Tg), which were identified using our novel oligo-capping-based protocol. Various types of cDNA data resources were interconnected with our original database functionalities. Specifically, in this update, we have included two unique RNA-Seq data sets consisting of 730 million mapped RNA-Seq tags. One is a dataset of 16 time-lapse experiments of cultured bradyzoite differentiation for Tg. The other dataset includes 31 clinical samples of Pf. Parasite RNA was extracted together with host human RNA, and the extracted mixed RNA was used for RNA sequencing, with the expectation that gene expression information from the host and parasite would be simultaneously represented. By providing the largest unique full-length cDNA and dynamic transcriptome data, Full-Parasites is useful for understanding host-parasite interactions and will help to eventually elucidate how monophyletic organisms have evolved to become parasites by adopting complex life cycles.

High-resolution characterization of Toxoplasma gondii transcriptome with a massive parallel sequencing method.

For the last couple of years, a method that permits the collection of precise positional information of transcriptional start sites (TSSs) together with digital information of the gene-expression levels in a high-throughput manner was established. We applied this novel method, 'tss-seq', to elucidate the transcriptome of tachyzoites of the Toxoplasma gondii, which resulted in the identification of 124,000 TSSs, and they were clustered into 10,000 transcription regions (TRs) with a statistics-based analysis. The TRs and annotated ORFs were paired, resulting in the identification of 30% of the TRs and 40% of the ORFs without their counterparts, which predicted undiscovered genes and stage-specific transcriptions, respectively. The massive data for TSSs make it possible to execute the first systematic analysis of the T. gondii core promoter structure, and the information showed that T. gondii utilized an initiator-like motif for their transcription in the major and novel motif, the downstream thymidine cluster, which was similar to the Y patch observed in plants. This encyclopaedic analysis also suggested that the TATA box, and the other well-known core promoter elements were hardly utilized.

DBTSS provides a tissue specific dynamic view of Transcription Start Sites.

DataBase of Transcription Start Sites (DBTSS) is a database which contains precise positional information for transcription start sites (TSSs) of eukaryotic mRNAs. In this update, we included 330 million new tags generated by massively sequencing the 5'-end of oligo-cap selected cDNAs in humans and mice. The tags were collected from normal fetal or adult human tissues, including brain, thymus, liver, kidney and heart, from 6 human cell lines in 21 diverse growth conditions as well as from mouse NIH3T3 cell line: altogether 31 different cell types or culture conditions are represented. This unprecedented increase in depth of data now allows DBTSS to faithfully represent the dynamically changing landscape of TSSs in different cell types and conditions, during development and in the course of evolution. Differential usage of alternative 5'-ends across cell types and conditions can be viewed in a series of new interfaces. Promoter sequence information is now displayed in a comparative genomics viewer where evolutionary turnover of the TSSs can be evaluated. DBTSS can be accessed at http://dbtss.hgc.jp/.

Evaluation of Echinococcus multilocularis tetraspanins as vaccine candidates against primary alveolar echinococcosis.

Echinococcus multilocularis causes an important zoonotic cestode disease. The metacestode stage proliferates in the liver of intermediate hosts including human and rodents and forms multiple cysts. Recently, members of a transmembrane protein tetraspanin (TSP) family have been used as vaccines against schistosomosis, or as diagnostic antigens for cysticercosis. In this study, seven tetraspanins of E. multilocularis, designated as TSP1 to TSP7, were evaluated for their protective potential against primary alveolar echinococcosis. The large extracellular loop (LEL) region of these tetraspanins was cloned from a full-length enriched cDNA library of E. multilocularis metacestodes and expressed in Escherichia coli as a fusion protein with thioredoxin. Recombinant TSPs were applied as vaccines against an E. multilocularis primary experimental infection in BALB/c mice. Cyst lesions in the livers of vaccinated and non-vaccinated mice were counted. The cyst lesion reduction rates induced by the seven tetraspanins in vaccinated vis-à-vis non-vaccinated mice were: 87.9%, 65.8%, 85.1%, 66.9%, 73.7%, 72.9% and 37.6%. Vaccination conferred protective rates to mice ranging from 0% (TSP5, 6, 7) to maximally 33% (TSP1, 3). The results indicated that recombinant tetraspanins have varying protective effects against primary alveolar echinococcosis and could be used in vaccine development.

Inconsistencies of genome annotations in apicomplexan parasites revealed by 5'-end-one-pass and full-length sequences of oligo-capped cDNAs.

BACKGROUND: Apicomplexan parasites are causative agents of various diseases including malaria and have been targets of extensive genomic sequencing. We generated 5'-EST collections for six apicomplexa parasites using our full-length oligo-capping cDNA library method. To improve upon the current genome annotations, as well as to validate the importance for physical cDNA clone resources, we generated a large-scale collection of full-length cDNAs for several apicomplexa parasites. RESULTS: In this study, we used a total of 61,056 5'-end-single-pass cDNA sequences from Plasmodium falciparum, P. vivax, P. yoelii, P. berghei, Cryptosporidium parvum, and Toxoplasma gondii. We compared these partially sequenced cDNA sequences with the currently annotated gene models and observed significant inconsistencies between the two datasets. In particular, we found that on average 14% of the exons in the current gene models were not supported by any cDNA evidence, and that 16% of the current gene models may contain at least one mis-annotation and should be re-evaluated. We also identified a large number of transcripts that had been previously unidentified. For 732 cDNAs in T. gondii, the entire sequences were determined in order to evaluate the annotated gene models at the complete full-length transcript level. We found that 41% of the T. gondii gene models contained at least one inconsistency. We also identified and confirmed by RT-PCR 140 previously unidentified transcripts found in the intergenic regions of the current gene annotations. We show that the majority of these discrepancies are due to questionable predictions of one or two extra exons in the upstream or downstream regions of the genes. CONCLUSION: Our data indicates that the current gene models are likely to still be incomplete and have much room for improvement. Our unique full-length cDNA information is especially useful for further refinement of the annotations for the genomes of apicomplexa parasites.